Multi-channel parallel/serial concatenated convolutional codes and trellis coded modulation encoder/decoder

ABSTRACT

A method and apparatus to transmit turbo-encoded data in a multitone channel assigns the original data and selected parity bits across multitone subchannels allowing transmission of an entire turbo-encoded block withing one or few symbol time frame. Parity bits are selected by a procedure using data derived by optimization using simulation of a single-channel system. The optimization determines, for a specified bit error rate, for each possible number of information bits per symbol, the code rate corresponding to the lowest signal-to-noise ratio. Alternatively, in the simulation non-identical integer values may be applied to the channels to approximate non-integer values of code rate and information bits-per-channel in the aggregate. The optimized data are used to determine an optimal code rate and SNR for each channel. In assigning the respective bits to the channels, the number of parity bits in each channel are shared.

The present application is a continuation-in-part of U. S. application Ser. No. 09/001,342 filed Dec. 31, 1997, now U.S. Pat. No. 6,088,387 entitled as above and by the same inventors.

FIELD OF THE INVENTION

The invention relates to turbo code (Parallel or Serial Concatenated Convolutional codes). More specifically, the invention relates to the combination of turbo code with trellis coded modulation (TCM) techniques within a multi-tone modulation environment.

BACKGROUND OF THE INVENTION

Several techniques for combining turbo code with trellis code modulation (TCM) are described in the current literature. The combination is attractive because it allows high spectral efficiency in the encoder/decoder to be combined with high coding gain. In all of the prior art approaches (see references 1 through 4), the main concern is achieving low bit error rate (BER) performance at the lowest signal-to-noise ratio (SNR) possible, while obtaining high spectral efficiency. The codes described in the prior art are designed for a fixed modulation system. For example, the codes are designed for M-Quadrature Amplitude modulation (M-QAM) or M-Phase Shift Keying (M-PSK), where M, the dimension of the constellation, is constant for the application.

Turbo code achieves very impressive bit error rate (BER) performance. However, turbo code exacts a price in terms of computational complexity and block size in exchange for such low BER. These problems are compounded when viewed from the standpoint of multichannel transmission systems. In multichannel or multicarrier systems, the number of bits that can be transmitted in a single symbol time in each channel is very small compared to the size of a turbo block. However, in order to decode turbo-encoded data, the entire block must be received at the modem. This is because all of the parity bits relate to all of the information bits. This means that the receiver must wait many symbol times in order to decode a single bit of data. For example, in QAM modulation in a multicarrier system, such as digital subscriber line (DSL) modems, where total power for all sub-bands is limited and a power spectral density (PSD) mask applies, QAM constellation sizes may be limited to something like 10 bits per channel. A typical turbo block size is of the order of 1000 bits plus the parity bits. The receiver would have to wait 150 symbol times (assuming a 2/3 code rate) and queue up a lot of data before it could decode.

In addition to the long delay for transmitting large turbo blocks, the prior art techniques, do not address the multi-tone modulation environments. In multi tone communications, multiple channels, each identified with a separate frequency carrier, are employed to transmit data. In a multi tone system, the channel characteristics vary from channel to channel and advantages are achieved by allowing the constellation sizes to vary from channel to channel; hence the name which is short for multiple tones. For example, determining how many parity bits to send is based on selecting a code rate that provides the highest possible information bit rate. The code rate, however, is a single parameter that must be traded against constellation size and signal-to-noise ratio (SNR). But in a multi tone environment, these are allowed to vary from channel to channel.

Since the filing of our prior U.S. patent application Ser. No. 09/001,342, filed Dec. 31, 1997, we have learned of two further prior art references, U.S. Pat. Nos. 5,446,747 and 5,563,897 which relate to Turbo code. The '747 patent describes two codings in parallel and the '897 patent is directed to detecting information bits processed by the codes.

SUMMARY OF THE INVENTION

An object of the invention is to provide an apparatus for implementing turbo code and trellis code modulation in a multi tone communications environment.

Another object of the invention is to provide an improved method for calculating turbo code parameters.

Still another object of the invention is to provide an improved method for coding and decoding turbo coded data bits.

Briefly, a method and apparatus to transmit turbo-encoded data in a multitone channel assigns the original data and selected parity bits across subchannels allowing transmission of an entire turbo-encoded block within one or few symbol time frames. Parity bits are selected by a procedure using data derived by optimization using simulation of a single-channel system. The optimization determines, for a specified bit error rate, for each possible number of information bits per symbol, the code rate and constellation size corresponding to the lowest signal-to-noise ratio. Alternatively, in the simulation non-identical integer values may be applied to the channels to approximate non-integer values of code rate and information bits-per-channel in the aggregate. The optimized data are used to determine an optimal code rate, constellation size, and SNR for each channel. In assigning the respective bits to the channels, the number of parity bits in each channel are shared

To eliminate the processing delay associated with transmitting large turbo blocks, this patent application proposes spreading a single turbo block across multiple channels. So, for example, if there are 512 channels in a multicarrier system and each channel carry 2 bits per symbol, the receiver can receive an entire turbo block plus parity bits in a couple of symbol times. In turbo code, the code rate and the constellation size can be varied selectively based on the channel characteristic. For example, turbo code produces at least two parity bits for each information bit. However, it is rare that the highest throughput is obtained by sending all the parity bits (the lowest code rate). The overhead is too high (that is, the code rate or ratio of information bits to total bits is too low) for the benefit of having so many parity bits. Thus, in using turbo codes, transmitters remove parity bits in a process of puncturing. Given a channel attenuation and noise level, one can find by simulation, the optimum number of parity bits to transmit and the optimum constellation size for the given power constraints. That is, constellation size is traded against signal-to-noise ratio and code rate for a fixed bit error rate and channel characteristic (attenuation and noise and whatever other characteristics may be folded into the channel model used for simulation) to optimize the number of information bits transmitted at the given BER.

The present invention applies a combination of turbo code and trellis code modulation (TCM) in a multi-tone modulation system. For example, the invention is applicable to a multi tone system employing M-quadrature amplitude modulation (M-QAM) or M-phase shift keying (M-PSK), where M, the dimension of the constellation, can vary from one channel to another channel. The turbo code trellis code modulation (TTCM) design of the invention, is based on the fact that the size of the constellation for each channel varies from channel to channel. The invention permits the use of a single TTCM encoder in the transmitter and a single TTCM decoder in the receiver. This feature reduces delay compared to prior art devices that employ separate individual encoders and decoders for each channel. The feature also reduces the complexity of the transmitter and receiver.

Another significant difference between the technique used in the invention and prior art techniques is the fact that bits are never estimated from the symbols. According to the invention, the received symbols are used to compute the parameters of the turbo code directly. In addition, in this regard, this technique is more efficient in terms of computational complexity. More specifically, combining log of the likelihood ratio approach and this technique makes computation of turbo code parameters very simple.

The invention can be used in systems where there are many available channels to transmit the information from one transmitter to one receiver. In such systems, instead of having multiple encoders and decoders at the transmitter and receiver to service all the channels, there will only be a single encoder that applies the turbo code trellis coded modulation to the channels and a single decoder to decode the information that is sent throughout the various channels. Any modulation scheme such as QAM or PSK can be used for this invention and the size of the constellation across all available channels can be different from each other. A similar technique can be used when we transmit the data over one channel, but the channel is divided into many frequency sub-bands and each one of these sub-bands is treated as a unique channel.

The current invention will only require one encoder for delivery of information over environment with multiple channels and decoding of the received encoded data with only one decoder at the receiver. The same advantages accrue where the channel is divided into many frequency sub-bands and each sub-band is treated as a separate channel.

BRIEF DESCRIPTION OF THE DRAWINGS

In the drawing,

FIG. 1 shows a transmitting portion of a multitone modem according to an embodiment of the invention having a channel assignment controller including a uniform interleaver.

FIG. 1A shows an alternate embodiment of the modem portion of FIG. 1 in which channel assignment controller 24 comprises parity bit selection and an external uniform interleaver is shown.

FIG. 1B shows a partial alternative embodiment, referred to herein as a self-concatenated code embodiment, of the modem portion of FIG. 1 in which RSC encoder 1 is combined with RSC encoder 2 and only a single RSC encoder remains providing a single parity bit stream to controller 24 via switch 15.

FIG. 2 shows a mapping scheme for a 16-QAM constellation.

FIG. 3 shows a receiving portion of a multitone modem according to an embodiment of the invention.

FIG. 3A shows a transmitting portion of a multitone modem according to another embodiment of the invention in which only the gamma values for the information bits are uniform-deinterleaved.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

Referring to FIG. 1, an encoder of a turbo-trellis code modulation (TTCM) system has first and second recursive systematic convolutional encoders (RSC encoders) 12 and 14. First RSC encoder 12 receives information bits from a data source 8 and generates a first stream of parity bits 17, one for each information bit from data source 8. Second RSC encoder 14 receives information bits from data source 8 through a random interleaver 16 and generates a second stream of parity bits 18, one for each information bit from data source 8. First and second RSC encoders 12 and 14 together with random interleaver 16 comprise a turbo encoder. Encoder types other than RSC encoders may be used in a turbo encoder, but, for purposes of the present invention, RSC encoders are preferred. In addition, these RSC encoders can be concatenated in series. Also, more than two RSC (or other type) encoders may be used in a turbo encoder, each receiving information through an interleaver. In the preferred embodiment, however, one such interleaver is used. Information bit stream 9 (data directly from data source 8) and parity bit streams 17 and 18 are applied to a channel assignment controller 24. In channel assignment controller 24, parity bits are selected from the parity bit streams 17 and 18 by a parity bit selection portion 21. The selected parity bits and information bit stream 9 are then interleaved by a second interleaver 22. FIG. 1 shows a single uniform interleaver 22 for information bits 9 and for receiving parity bits selected by parity bit selection 21, FIG. 1A shows uniform interleaver 22 for uniformly interleaving information bits 9 and providing the interleaved information bits to controller 24. The second interleaver 22 can actually include two interleavers. If two interleavers are used, one interleaver is for information bits and the second one is for parity bits. Second interleavers can be either random or uniform interleavers. In this application uniform interleavers are preferred. The interleaved bits from uniform interleaver 22 are applied to a QAM modulator 26. Channel assignment controller 24 sends signal power data 27A to QAM modulator 26 to control the power level of the respective QAM tones generated by it. Also, Channel assignment controller 24 sends mapping data 27B to QAM modulator 26 to assign each bit in the interleaved data stream to respective coordinates in the respective QAM tones generated by QAM modulator 26. QAM modulator outputs QAM tones to a channel. The details of how channel assignment controller 24 performs the tasks identified above are discussed below. Note that QAM modulator 26 could be replaced by a phase shift keying (PSK) modulator or any other modulation method could also be used.

Channel assignment controller 24 performs the selection and direction of information and parity bits according a specific set of rules that uniquely determines the constellation size for each channel and which information and parity bits are assigned to each QAM coordinate. In addition, channel assignment controller 24 also directs the power level for each QAM tone generated in each channel by QAM modulator 26. These features of channel assignment controller 24 are described below.

Referring to FIG. 1B, a simplified embodiment portion which may be referred to as a self-concatenated code embodiment of FIG. 1 is shown in which a single RSC encoder 14 replaces the two RSC encoders 1 and 2 shown. In this embodiment, a data source 8 provides information bits 9 to a random interleaver 16 and to RSC encoder 14 via a switch 15. Periodically, switch 15 switches from information bits 9 to interleaved information bits and RSC encoder provides a parity bit stream in time sequence from both switched sources. Controller 24 receives a first parity bit stream portion and stores that portion until the second parity stream portion is received. Thus, controller 24 is then able to select parity bits from each stream as in a two RSC encoder embodiment. The single RSC encoder 14 may be said to be time shared to produce both parity bit streams as one stream 18 in time sequence. Random interleaver 16 provides an interleaved information bit stream to RSC encoder 14 via switch 15 and data source 8 provides information bit stream 9 via switch 15 in time sequence according to which switch position switch 15 is in. The single RSC encoder has the reference numeral 14 for consistency with FIG. 1 and “RSC encoder 1” 12 is eliminated. Because encoder 12 is eliminated, RSC encoder 14 provides only a single time-shared parity bit stream 18 to channel assignment controller 24. Channel assignment controller 24 may have a uniform interleaver 22 contained therein or external as shown in FIG. 1 or 1A. Parity bit selection 21 of controller 24 operates on the single time-shared parity stream by selecting from one parity bit stream sequence stored in controller memory (not shown) or the other switch parity bit stream sequence as if it received two separate streams from two separate RSC encoders to select parity bits for transmission.

Referring to FIG. 1, QAM tones generated by QAM modulator 26 are applied to a single physical channel with multiple channels. The respective QAM tones are received by a QAM demodulator 30 shown in FIG. 3. QAM demodulator 30 separates the I AND Q channel data for each QAM tone. These are applied to a gamma function processor 32 which calculates a gamma function from the I and Q channel data. Gamma function data represents a first order probability of both the information bits and the parity bits and its derived based on the channel characteristics and mapping of the information bits and parity bits into a QAM signal. The gamma probabilities are calculated in real time as I and Q channel data are sequentially generated from the multi tone signal. The gamma probabilities are applied to a uniform deinterleaver 28 which reverses the interleaving performed by uniformed interleaver 22. The data resulting from the latter deinterleaving process represents, probabilistically, the original (punctured) parity bit stream and the original information bit stream in the order in which it was delivered from the turbo-encoder prior to uniform interleaving. For the gamma functions to be calculated properly in the first instance by gamma function processor 32, gamma function processor 32 must receive mapping data 48 which gives constellation size and bit assignments (i.e, which bits are coded by the respective QAM tones for each channel). This information is provided by channel assignment controller 24 in the receiving modem. Since information is provided derived during initialization is provided to both modems and since the assignment of data from turbo encoder 10 is uniquely determined by rules programmed into channel assignment controller 24, channel assignment controller 24 is capable of providing control information to gamma function processor 32 required to calculate correct gamma probabilities from the I and Q channel data. The deinterleaved data are fed to a turbo decoder 49 for iterative decoding using an appropriate turbo-decoding process, such as described in the literature.

The number of bits (m) transmitted in each channel depends on the power of noise, the channel attenuation, PSD mask limit and any other characteristics of that particular channel that may vary with frequency. The number of bits, including both information and parity bits, corresponds to the constellation size of an M-QAM or M-PSK symbol. If there are m bits to be represented by a QAM or PSK tone, there must be M=2^(m) different distinguishable QAM or PSK symbols and the constellation size is said to be M. During the initialization procedure a pseudo-random signal is transmitted at one end of the channel and received at the other end to measure the channel transform characteristic. The transform characteristic is used to determine the respective constellation size for each channel using the channel data stored as a table showing the channels and corresponding parameter(s), for example, the noise power-gain ratio, for each channel.

The size of random interleaver 16 corresponds to the size of the block that is chosen for the turbo code in the particular application. Uniform interleaver 22 can be smaller than the first interleaver. The size of uniform interleaver 22 should consider factors such as the duration of the impulsive noise in the channel, or, if there is a fading channel, its size is based on the maximum duration of the fade in all the channels. In general, it is possible to interleave only the information bits in uniform interleaver 22 without interleaving the parity bits together with the information bits as indicated in FIGS. 1A and 3A. As indicated in FIG. 3A, the parity bits are not uniformly deinterleaved because this is not necessary since the parity bits are not uniformly interleaved by the transmitter of FIG. 1A.

Referring to FIG. 2, the assignment of bits to the (e.g., M-QAM) symbols can employ any suitable mapping scheme such as Gray mapping, natural mapping, or reordered mapping. One of the rules by which channel assignment controller 24 controls the mapping of bit sequences to symbol coordinates involves the concept of a most-protected bit. In FIG. 2, a 16-QAM signal mapping is shown. As mentioned, a 16-QAM signal can represent 4 bits and each coordinate shown in the figure represents a possible value for a 4-bit word represented by the QAM signal. Of the four bits 41 through 44, the two most significant bits 41 and 42 are most protected. As can be seen in FIG. 2, most significant bit 41 has the value of one for all points in an envelope 46 (and consequently it has the value of zero for all the other points). Similarly, the second most significant bit 42 has the value of one for all points in envelope 45 (and consequently it has the value of zero for all the other points). These bits are the most protected because the likelihood is that the largest degree of distortion in the QAM tone is required for one of these bits to be incorrectly represented. Note that for second most significant bit 42 to be misrepresented, there must be distortion in the I channel signal. For most significant 41 to be misrepresented, there must be distortion in the Q channel. For many mapping schemes, some bits will have this property of being less likely for a given level of distortion to be misrepresented in the final distorted signal. For QAM and PSK and other kinds of symbols, there will be two most protected bits, one corresponding to the I channel and one corresponding to the Q channel.

In assigning bits to the various symbols, channel assignment controller 24 adheres to a number of rules. A first rule according to which channel assignment controller 24 assigns bits governs the selection of parity bits to preserve and transmit (the others being rejected). As mentioned above, the number of parity bits created by the turbo-encoder is far more than is justified to transmit. The benefit of transmitting a large number of bits for error-correction is offset by the consequent reduction in the total number of information bits since the communication system is only capable, to an approximation, transmitting a certain total number of bits (information bits plus parity bits) per unit time for a given total transmit power. The optimal number of parity bits depends on the channel characteristics and limits on power. Optimal code rates of more than ½ and higher are typical. Therefore, in most cases, many parity bits are simply not transmitted. The number of parity bits transmitted should therefore be optimized based on the channel characteristics.

In a network comprising multiple identical channels, the code rate can be optimized by simulation. The goal of the optimization is to determine, for a given BER and a given number of information bits, the combination of the total number of bits and the total number of parity bits that corresponds to the lowest total SNR An AWGN (or other type of) channel can be simulated numerically (or a real-world simulation using a pseudo-noise signal in a real channel could be done). Since there is a single constellation size in this situation, the code rate can be traded against signal power and constellation size to arrive at an optimum code rate for a predefined BER. In the multi tone environment, the constellation size may vary across many channels. To simulate every possible spectrum of channel characteristics would be impractical. Moreover, the amount of memory required to store results of such an optimization, to permit implementation, would also be impractical and undesirably expensive.

It has been found through simulation that optimization results from a single-channel network can be extended to channels with variable channel characteristics according to a procedure which is now described. Near-optimum performance in the multi tone environment can be obtained by optimizing the constellation size, code rate (i.e., number of parity bits), and SNR for a single channel whose bandwidth corresponds to the bandwidth of the channels in a multi tone network. A simulation model that calculates SNR given a required minimum BER, constellation size, number of parity bits (i.e., code rate) is used to determine, for a given BER and number of information bits, the combination of total number of bits and parity bits that correspond to the lowest SNR. The result will be a table which specifies, for each value of information bits per channel, an optimal value of SNR, constellation size, and parity bits per channel (the latter could be specified as code rate).

Note that to perform the above optimization, a driving program, such as available from math libraries, can be used to alter the various parameters in a way that locates the optimized points with as few simulations as possible. For each number of information bits per channel, an optimal combination of constellation size, SNR, and number of parity bits is found. The optimization data is stored in a table (called the optimized rate table or ORT) and used by channel assignment controller 24 to direct the assignment of bits to the various symbol coordinates for each channel. The results of the single-channel optimization are used to optimize the code rate in the multi tone environment as described below.

In the ORT, each row corresponds to one of the possible number of information bits and the corresponding optimum SNR and code rate (the constellation size being determined by the number of information bits and the code rate). Channel assignment controller 24 assigns bits to each channel in turn using the ORT data and the channel characteristic data measured during the initialization procedure. Channels are preferably sorted by a channel characteristic parameter, such as the difference between PSD mask limit and the ratio of gain to noise power, and assigned bits sequentially through the resulting series. The first step in assigning bits for a particular channel is to determine the total number of bits to be transmitted in the channel and the number of those bits that will be parity bits. This is determined by first assigning a constellation size and number of parity bits for each channel. Working in descending order by number of information bits the SNR data in the ORT is compared to the channel characteristic and any PSD mask limit that applies to the given channel. If the required power exceeds the limit, the next lower number of information bits is checked until the total power required falls below the mask limit. This provides a constellation size, code rate, and SNR for the channel.

An alternative procedure for generating the ORT and implementing channel assignment may be used to obtain a finer degree of resolution in the independent parameters. Consider what happens when a simulation determining the optimum number of parity bits per channel is based on a single channel (or viewed as an optimum for each channel alone). This forces one to select the total number of parity bits to be transmitted, when speaking of multiple channels, from integral multiples of the number of channels. Suppose there are 256 channels and the turbo block as 1024 bits in it. Suppose further that 128 of the channels have very similar transfer characteristics. If a single channel model is used to determine the optimum for these 128 channels, the optimum number of parity bits for the entire turbo block as a whole must be selected from among multiples of 128 bits. The modem can send 1 parity bit per channel or 2 parity bits per channel, etc. For these 128 channels, the code rate is varying by huge degrees with each change in the number of parity bits. It must be varied in chunks of 128 bits. For just this group of channels, that corresponds to a step size in code rate of 0.125. The number of information or parity bits per channel may be permitted to vary from channel to channel, even though they have nearly identical transfer characteristic (identical for the simulation and nearly identical for the practical situation). Thus, consider a group of five channels and ask, for a given number of total information bits in these five channels, what is the total number of parity bits (or code rate) and constellation size that corresponds to the lowest signal to noise ratio? The optimum might pack a total of 8 parity bits in the five channels, which corresponds to some channels carrying a parity bit and some channels carrying no parity bit. For example, consider 3.33 information bits per channel and a total of 4 bits per channel. The bits could be loaded for the simulation as follows:

Infor- Channel mation Parity Total 1 3 1 4 2 3 1 4 3 3 1 4 4 4 0 4 5 3 1 4 6 3 1 4 7 4 0 4 8 3 1 4 9 3 1 4 10 4 0 4 11 3 1 4 12 3 1 4 13 3 1 4 14 4 0 4 15 3 1 4 16 3 1 4 17 4 0 4 18 3 1 4

Since each channel can carry only an integral number of parity bits and information bits, to vary total code rate (summed over all channels) the number of parity bits must be varied from channel to channel. A scheduling formula to allocate parity bits for purposes of simulation among the channels may take the form (although other scheduling formulas are possible): $\begin{matrix} {p_{i} = {{{INT}\left\lbrack {{\sum\limits_{j = 1}^{i}\quad P_{j}} - {\sum\limits_{k = 1}^{i - 1}\quad p_{k}}} \right\rbrack} = {{INT}\left\lbrack {\left( {{\sum\limits_{j = 1}^{i - 1}\quad P_{j}} - p_{j}} \right) + P_{i}} \right\rbrack}}} & (1) \end{matrix}$

where p_(i) is the number of parity bits in the ith channel, and P_(i) is the non-integer value of the number of parity bits per channel for the ith channel. For purposes of simulation, the values P_(i) are the same for all channels (The same equation is used to assign integer values to the multi tone channels and in that case, the value of P_(i) varies from channel to channel). The above equation will produce a series of integers that approach, on average, the non-integral values for a large number of channels.

The following procedure may be used to perform the simulation. Instead of simulating a single channel and optimizing the single channel, a number of identical channels are simulated, that is, they all have the same frequency characteristic. This optimization would produce exactly the same result as a single channel if each channel were required to carry the exact same number of information bits and parity bits. However, in this case, the number of parity bits and information bits are allowed to vary from channel to channel. So, considering a set of channels, the optimization does not seek to determine, for a specified bit error rate, for each possible number of information bits per symbol, the code rate and constellation size corresponding to the lowest signal-to-noise ratio. Instead, the optimization looks at the channels as a group and asks a similar question. It could be, for a given number of information bits per channel (not necessarily integer) what is the constellation size that corresponds to the lowest signal to noise ratio? That is, the channels are viewed as a group, rather than optimized individually. Thus, for example, to simulate 3.5 information bits per channel, the constellation size may be stepped through a range of sizes from 4 to, say, 8, the constellation size being the same for all the channels for each simulation run. Every other channel would have three information bits and the rest would have 4 information bits. For a constellation size of 4 the number of parity bits would alternate between 0 and 1 in adjacent channels so that the average number was 4-3.5 or 0.5 parity bits per channel. The rest of the constellation sizes would be simulated and the result would be a table as discussed above, except that non-integer values of parity bits per channel or information bits per channel would be allowed.

Equation (1) above applies to the situation where the channel characteristic of all frequency channels is identical. To apply non-integer code rates and/or information bits per channel to each channel in a process similar to that defined by equation (1) except that the values of the value of P_(i) varies from channel to channel. To make the assignment of bits to the channels approach the situation simulated, where adjacent channels had identical characteristics, in performing the bit allocation, the frequency bins may be sorted by some channel characteristic so that the optimal parameters vary slowly and smoothly as bits are assigned in sequence. Thus, for example, in assigning bits, an adjacent series of channels would tend to be similar enough to have the same optimal number of information bits per channel allocated to them. As the bits are assigned, a sequence equal to the non-integer loading factor would be maintained in applying equation (1). If the number of bits per channel is allowed to vary abruptly, the assignment of bits would not relate well to the simulation that determined the optimal values.

Once parity and information bits are assigned to all the bins, if the aggregate power required to transmit all the bits in all the channels exceeds any applicable frequency-aggregate power limit, bits can be removed from the bins. Preferably, to meet the aggregate power limit, information bits are removed from the bins by determining the step in the ORT toward lower SNR corresponding to the greatest reduction in signal power in exchange for the least amount of information change. This can be done using the ORT in a procedure described in U.S. Patent Application entitled Method and Apparatus for Allocating Data for Transmission Via Discrete Multiple Tones- Sonalkar, et. al, the entirety of which is incorporated herein by reference. The only modification to the procedure described in the above application is that, instead of removing bits per se, the procedure uses the ORT table to determine the power required to transmit data rather than a formula. In this modified procedure, the algorithm would step through the column of SNRs to fill up the bins to the PSD mask limit. Then, to satisfy the aggregate power limit, if necessary, the procedure would go through the ORT to determine the incremental change in constellation/SNR assignment for a channel that would reduce the aggregate power the most for the smallest change in total information bits transmitted. In essence, the energy calculation function described in that reference is replaced by the SNR determined by the ORT.

Although in the preferred embodiment described above, the optimized rate table (ORT) is derived using a single-channel model, other possible optimization methods are possible. An alternative optimization method which can be used to generate the ORT data is to optimize the total number of information bits transmitted for a series of representative profiles of channel characteristic with frequency. These results would provide an optimal code rate for the entire series of channels. A best fit between the representative frequency profiles and the measured profile obtained during initialization would indicate the code rate for the transmission. Then, using the code rate, each channel could contribute to the transmission of parity bits on a pro rata basis according to constellation size determined by a conventional bit-loading algorithm used in multi tone systems.

Exactly which parity bits are preserved by transmission with the information bits and which are rejected is determined by another set of rules. In addition, the assignment of the respective bits to the respective symbol coordinates is based on still other rules. These rules are described below.

To determine which of the parity bits to allocate for transmission, channel assignment controller 24 selects parity bits originating from among the multiple RSCs such that, on average, each RSC encoder will contribute an equal the number of parity bits to the total transmitted over all channels for a given block. In addition, the parity bits are assigned in such a way as to distribute parity bits from each RSC encoder as uniformly across the spectrum as possible. This is done by allocating parity bits sequentially and alternating among the RSCs.

Another rule followed by channel assignment controller 24 relates to the assignment of information and parity bits to specific symbol coordinates. If there is one parity bit to be represented in a channel, the parity bit is assigned to one of the most protected bits. To insure that the assignment is unique, the I-channel most protected bit can be preferred or the controller can alternate between the I-channel most protected bit and the Q-channel most protected bit. If there are two parity bits to be conveyed in a channel, the parity bits are split between the I-channel most protected bit and the Q-channel most protected bit. If there are three parity bits to be assigned, the first and second most protected bits will be assigned. It would be rare for there to be more than 2 parity bits in a 16-QAM constellation as code rates for turbo-coded signals are usually higher than ½. Where no n^(th)-most protected bit can be identified, some rule that makes the choice for assigning any remaining parity bits can be chosen as long as it is unambiguous and uniquely determines the assignment of the parity bits.

When the end of a turbo block is reached, the controller treats any remaining space available in the current bin as a new channel and continues the allocation procedure beginning with the first channel in the next symbol period (as if it were wrapped around to the current symbol period). When the last channel is allocated in a symbol period, the symbols are generated.

Receiver 2 and the transmitter 1 both either have the information about the symbol constellation sizes a priori or they share the information during a startup or refresh period. Since both the receiver and the transmitter both have the same information, by applying the procedures and rules discussed above, their respective channel assignment controllers 24 can independently determine which bits are assigned to the various symbol coordinates.

The invention employs a unique decoding process. To explain this in detail, it is necessary to define some of the parameters of the turbo code. Each symbol (for example if a QAM symbol is used) may be represented as a complex number at the output of TTCM encoder 10. This complex number, A_(k)+jB_(k),is transmitted through the channel and at the receiver, additive noise in the channel is combined with the original signal to produce a new signal represented by the following equations for the I and Q channels:

I _(k) =A _(k) +n _(i)  (2)

Q _(k) =B _(k) +n _(q)  (3)

where n_(j), and n_(q) are additive noise components, in general, not necessarily Gaussian noise. In prior art systems, the symbols are first soft-decision decoded by computing the Logarithm of Likelihood Ratio (LLR). The soft-decoded data are applied to an input to a turbo decoder. The LLR is computed as follows: $\begin{matrix} {{{\lambda \quad \left( {x\quad (k)} \right)} = {K\quad {{Log}\left\lbrack \frac{\left. {{{\Pr \left\{ {x\quad (k)} \right)} = {1/I_{k}}},Q_{k}} \right\}}{\left. {{{\Pr \left\{ {x\quad (k)} \right)} = {0/I_{k}}},Q_{k}} \right\}} \right\rbrack}}},{i = 1},\ldots \quad,{m - m^{\prime}}} & (4) \end{matrix}$

where K is a constant, x(k) is the actual value of the information bit transmitted. Decoding of parity bits is performed by the same equation, replacing x(k) in (4) by c_(i)(k) where I=1,2, . . . , L−1, where (L−1) is the total number of RSC encoders. In some papers, a simplified version of equation (4) is suggested. In our approach, we utilize the I and Q channel data, I_(k) and Q_(k), directly to calculate the parameters of the turbo code. Before explaining this procedure, it is necessary to describe the decoding procedure in the turbo code briefly.

Assume that the turbo encoder is based on the parallel concatenation of RSC encoders (the RSC encoders can also be concatenated in series) and that the length of the shift register in each RSC encoder is n. In general, the number of shift registers for each RSC encoder can be different. For the case where the number of shift registers, n, is the same for all of the RSC encoders, there will be 2^(n)=N different states for each one of the RSC encoders. The LLR estimate of each decoded bit can be written based on the a-posteriori probability (APP) of that bit. $\begin{matrix} {{\lambda \quad \left( {x\quad (k)} \right)} = {K\quad {{Log}\left\lbrack \frac{\left. {{\Pr \left\{ {x\quad (k)} \right)} = {1/{observation}}} \right\}}{\left. {{\Pr \left\{ {x\quad (k)} \right)} = {0/{observation}}} \right\}} \right\rbrack}}} & (5) \end{matrix}$

The state of encoder 10 at time k is defined as S_(k) and the sequence of data (observation) is defined as R₁ ^(K)={R₁, . . . , R_(K)}. Each observation datum is defined as R_(k)={x(k), ci(k)} where x(k) and c_(i)(k) represent the received information bit and the parity bit of the I-th encoder at time k respectively after noise is added to the signal. The parity bit may not be available for some of the observation data. The number of RSC encoders is assumed to be 2, in this example. The block length of the data sequence is K bits. The APP of a decoded bit x(k) can be written as:

Pr {x(k)=i/R ₁ ^(K}=ΣΨ) _(k) ^(i) (N), i=0, 1  (6)

where Ψ_(k) ^(i)(N)·is defined as:

Ψ_(k) ^(i) (N)=Pr {x(k)=i, S _(k) =N/R ₁ ^(K)}.  (7)

From (6) and (7), (5) can be written as follows: $\begin{matrix} {\lambda_{1} = {{Log}\quad \frac{\sum\limits_{N}^{\quad}\quad {\psi_{k}^{1}\quad (N)}}{\sum\limits_{N}^{\quad}\quad {\psi_{k}^{0}\quad (N)}}}} & (8) \end{matrix}$

Ψ_(k) ^(I)(N) in (8) can be written in terms of three parameters. These parameters are defined as: $\begin{matrix} {{\alpha_{k}\quad (N)} = {{\Pr \left\{ {S_{k} = {N/R_{1}^{k}}} \right\}} = {h_{\alpha}\quad {\sum\limits_{N1}^{\quad}\quad {\sum\limits_{i = 0}^{1}\quad {\gamma_{i}\quad \left( {R_{k},{N1},N} \right)\quad \alpha_{k - 1}\quad ({N1})}}}}}} & (9) \\ {{\beta_{k}\quad (N)} = {\frac{\Pr \left\{ {{R_{k + 1}^{K}/S_{k}} = N} \right\}}{\Pr \left\{ {R_{k + 1}^{K}/R_{1}^{K}} \right\}} = {h_{\beta}\quad {\sum\limits_{N1}^{\quad}\quad {\sum\limits_{i = 0}^{1}\quad {\gamma_{i}\quad \left( {R_{k + 1},{N1},N} \right)\quad \beta_{k - 1}\quad ({N1})}}}}}} & (10) \end{matrix}$

 γ_(i)(R _(k) , N 1, N)=Pr{x(k)=i, S _(k) =N, R _(k) /S _(k−1) =N}  (11)

γ_(i)(R _(k) , N 1, N)=Pr{R _(k) /x(k)=i, S _(k) =N, S _(k−1) =N 1 }·Pr{x(k)=i/S _(k) =N, S _(k−1) =N 1}·Pr{S _(k) =N/S _(k−1)=1}  (12)

where h_(α) and h_(β) are the normalization factors. Ψ_(k) ^(i)(N) can then be written as: $\begin{matrix} {{{\psi_{k}^{i}\quad (N)} = {\sum\limits_{N1}^{\quad}\quad {\gamma_{i}\quad {\left( {R_{k},{N1},N} \right) \cdot \alpha_{k - 1}}\quad {({N1}) \cdot \beta_{k}}\quad (N)}}},{i = 0},1} & (13) \end{matrix}$

Then (8) can be re-written using (13). $\begin{matrix} {{\gamma_{1}\quad \left( {x\quad (k)} \right)} = {{Log}\quad \frac{\sum\limits_{N}^{\quad}\quad {\sum\limits_{N1}^{\quad}\quad {\gamma_{i}\quad {\left( {R_{k},{N1},N} \right) \cdot \alpha_{k - 1}}\quad {({N1}) \cdot \beta_{k}}\quad (N)}}}{\sum\limits_{N}^{\quad}\quad {\sum\limits_{N1}^{\quad}\quad {\gamma_{0}\quad {\left( {R_{k},{N1},N} \right) \cdot \alpha_{k - 1}}\quad {({N1}) \cdot \beta_{k}}\quad (N)}}}}} & (14) \end{matrix}$

In an Additive White Gaussian Noise (AWGN) channel, the first factor in the right hand side of equation (12), Pr{R_(k)/ . . . }, is the transition probability of the discrete Gaussian memoryless channel. Equation (12) can be determined from transition probabilities of the discrete memoryless channel and transition probabilities of the encoder trellis. Assuming the condition (x(k)=I, S_(k)=N, S_(k−1)=NI), R_(k) consists of two uncorrelated Gaussian variables, namely x(k) and ci(k), and thus we can write it as the product of two probabilities:

Pr{R _(k) /x(k)=i, S _(k) =N, S _(k−1) =N}=Pr{x _(r)(k)/ . . . }·Pr{ci _(r,i)(k)/ . . . }  (15)

In many applications, the x(k) and ci(k) are softly estimated by calculating LLR using equation (4) or a simplified version of (4). Then, the result will be used to find the conditional probabilities defined in (15). However, since the receiver is aware of the constellation of the signals, the symbols can be used directly to compute the probabilities in (15).

Pr{x _(r)(k)/x(k)=i, S _(k) =N, S _(k−1) =N1}=Pr{Ik, Qk/ . . . }

Pr{c _(r,i)(k)/x(k)=i, S _(k) =N, S _(k−1) =N1}=Pr{Ik, Qk/ . . . }  (16)

Equation (16) is used to compute the gamma function (the first parameter in the right hand side of (12)). This approach reduces the computational complexity of the system without degrading the performance of the system. After computing the gamma function, the alpha and beta functions can be computed according to equations (10) and (11). If this approach is combined with logarithm of the likelihood ratio approach in which the logaritm of the turbo parameters are calculated, then for an AWGN channel calculation of the gamma function will be substantially simplified.

The respective entireties of the following publications are incorporated herein by reference.

[1] Stephane Le Goff, Alain Glvieux, and Claude Berrou, “Turbo-Codes and High Spectral Efficiency modulation” Proceedings of ICC'94, New Orleans, La., May 1994.

[2] Patrick Robertson and Thomas Worz, “A novel bandwidth efficient coding scheme employing turbo-codes”, Proceedings of ICC'96, Dallas, Tex., June 1996.

[3] Sergio Benedetto, Dariush Divsalar, Guido Montorsi, and Fabrizio Pollara, “Parallel concatenated trellis coded modulation”, Proceedings of ICC'96, Dallas, Tex., June 1996.

[4] Benedetto, Divsalar, Montorsi, Pollara, “Serial Concatenated Trellis Coded modulation with Iterative Decoding”, ISIT97.

[5] Claude Berrou, Alain Glvieux, and Punya Thitimajshima, “Near Shannon Limit Error Correcting Coding and Decoding: turbo-codes”, Proceedings of ICC'93, Geneva, Switzerland, May 1993.

[6] Claude Berrou, and Alain Glvieux, “Near Optimum Error Correcting Coding and Decoding: turbo-codes,” IEEE Transaction on Communications, October 1996.

Our prior similarly titled patent application, U.S. patent application Ser. No. 09/001,342, filed Dec. 31, 1997, should be deemed to be incorporated by reference for all purposes including patent applications identified therein to be incorporated by reference as to their entire contents. 

We claim:
 1. A data encoder circuit for use in a data modem for communicating over a network having channels comprising: a random interleaver for receiving information bits from a data source and outputting interleaved information bits; a switch for receiving information bits from said data source and said output interleaved information bits and for switching between said information bits and said output interleaved information bits and outputting a resultant information bit stream; a recursive systematic convolutional encoder receiving the resultant information bit stream and outputting a parity bit stream comprising first and second time sequenced portions; and a controller comprising a parity bit selector and a memory, the controller storing said first time sequenced portion of said output parity bit stream from said recursive convolutional encoder in memory and, upon reception of said second time sequenced portion of said output parity bit selecting party bits for transmission.
 2. A method of encoding data for use in a data modem for communicating over a network having channels comprising the steps of: randomly interleaving information bits received from a data source and outputting interleaved information bits; periodically switching between information bits from said data source and said output interleaved information bits and outputting a resultant information bit stream; encoding via a recursive systematic convolutional encoder, the encoder receiving the resultant information bit stream and outputting a parity bit stream comprising first and second storing said first time sequenced portions; storing said first time sequenced portion of said output parity bit stream from said recursive convolutional encoder in memory; and upon reception of said second time sequenced portion of said output parity bit stream, selecting parity bits for transmission. 